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Annihilation of superheavy particles in primordial bound systems is considered as the source of Ultra High 
Energy Cosmic Rays (UHECR). Charge conservation makes them to be produced in pairs, and the estimated 
separation of particle and antiparticle in such pair is shown to be in some cases much smaller than the average 
separation determined by the averaged number density of considered particles. If the new U(l) charge is the source 
of a long range field similar to electromagnetic field, the particle and antiparticle, possessing that charge, can 
form primordial bound system with annihilation timescale, which can satisfy the conditions, assumed for this type 
of UHECR sources. These conditions severely constrain the possible properties of considered particles. So, the 
proposed mechanism of UHECR origin is impossible to realise, if the U(l) charged particles share ordinary weak, 
strong or electromagnetic interactions. It makes the proposed mechanism of pairing and binding of superheavy 
U(l) charged particles an effective theoretical tool in the probes of the physics of very early Universe and of the 
hidden sector of particle theory, underlying it. 



1. INTRODUCTION 

Ultra High Energy Cosmic Rays (UHECR) is 
the observed effect of superhigh energy physics 
in the modern Universe, what naturally puts to- 
gether particle physics and cosmology in the anal- 
ysis of their possible origin and effects. Particle 
theory predicts a wide variety of new phenom- 
ena in this energy range. Such phenomena are 
unavoidable in the modern Big bang cosmology, 
based on inflationary models with baryosynthcsis 
and (multicomponent?) nonbaryonic dark mat- 
ter. The physics of inflation and baryosynthesis, 
as well as dark matter and energy content implies 
new particles, fields and mechanisms, predicted in 



the hidden sector of particle theory. Such parti- 
cles, fields and mechanisms may play an impor- 
tant role in the problem of UHECR. It makes 
new physics necessary component of the analysis 
of UHECR data. 

The origin of cosmic rays with energies, exceed- 
ing the GZK cut off energy pQ, is widely discussed 
as the possible effect of new physics. One of pop- 
ular approaches is related with decays or anni- 
hilation in the Galaxy of primordial super-heavy 
particles |2l:ij (see 0] for review and references 
where in). 

The mass of such super-heavy particles to be 
considered in present paper is assumed to be 
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higher than the re-heating temperature of infla- 
tionary Universe, so it is assumed that the parti- 
cles are created in some non-equilibrium processes 
(see e.g. |5I6| 1. taking place after inflation at the 
stage of preheating. 

The problems, related with this approach, are 
as follows. If the source of ultra high energy cos- 
mic rays (UHECR) is related with particle decay 
in the Galaxy, the timescale of this decay, which is 
necessary to reproduce the UHECR data, needs 
special nontrivial explanation. Indeed, the relic 
unstable particle should survive to the present 
time, and having the mass m of the order of 10 14 
GeV or larger it should have the lifetime r, ex- 
ceeding the age of the Universe. On the other 
hand, even, if particle decay is due to gravita- 
tional interaction, and its probability is of the or- 
der of (here and further, if not directly indicated 
otherwise, we use the units h = c = k = 1) 



where Mp = 10 19 GeV is the Planck mass, the 
estimated lifetime would be by many orders of 
magnitude smaller. It implies some specific ad- 
ditional suppression factors in the probability of 
decay, which need rather nontrivial physical real- 
ization ( |2I4| ). e.g. in the model of cryptons |7] 
(see jS] for review). 

If the considered particles are absolutely sta- 
ble, the source of UHECRs is related with their 
annihilation in the Galaxy. But their averaged 
number density, constrained by the upper limit 
on their total density, is so low, that strongly in- 
homogeneous distribution is needed to enhance 
the effect of annihilation to the level, desired to 
explain the origin of UHECR by this mechanism. 

In the present paper, we consider in more de- 
tails new approach to the solution of the latter 
problem, offered in [JJJ. If superheavy particles 
possess new U(l) gauge charge, related to the hid- 
den sector of particle theory, they are created in 
pairs. The Coulomb-like attraction (mediated by 
the massless U(l) gauge boson) between parti- 
cles and antiparticles in these pairs can lead to 
their primordial binding, so that the annihilation 
in the bound system provides the mechanism for 
UHECR origin. To realize this mechanism the 



properties of superheavy particles should satisfy 
a set of conditions, putting severe constrains on 
the cosmological scenarios and particle models, 
underlying the proposed mechanism. 

The necessary decoupling of superheavy par- 
ticles from the interactions of ordinary particles 
can be related with physics of neutrino mass, re- 
sulting in the dominant annihilation channels to 
neutrino. It may be important for another ap- 
proach to a possible solution of the GZK para- 
dox that considers the Ultra High Energy Cos- 
mic Rays as secondary products of UHE neu- 
trinos, originated at far cosmic distances, over- 
coming GZK cut-off, hitting onto relic light neu- 
trino in Hot Dark Halos, leading to resonant Z 
boson production. A consequent Z-Shower (Z- 
Burst)(see |l()lllll2li:flT3) ) takes place, where a 
boosted ultra-relativistic gauge boson Z (or WW, 
ZZ pairs) decays in flight and where its UHE 
nuclear secondaries are the observed UHECR 
events in terrestrial atmosphere. These ZeV pri- 
mary UHE neutrinos may be produced either in- 
side compact astrophysical objects (Jets GRBs, 
AGNs,BL Lac 15_ ) or by relic topological defects 
decay or, as in the present paper, by super- 
heavy particle annihilation in primordial bound 
systems. In the first (compact object) case one 
may easily understand the observed UHECR clus- 
tering as well as the possible correlation found re- 
cently with BL Lac sources. In the second case 
one may assume the UHECR clustering toward 
BL Lac as a pure coincidence; otherwise, one may 
consider a possible faster induced annihilation of 
these bound systems inside deeper clustered grav- 
itational wells around AGN and BL Lacs objects, 
treating such annihilation as the source of UHE 
neutrinos. 

2. Non-equilibrium mechanisms of super- 
heavy particles production 

The modern cosmological models of inflation- 
ary Universe assume that thermodynamically 
equilibrium conditions of hot Universe (the so 
called "reheating") do not take place immedi- 
ately after the end of inflation, and that there ex- 
ist rather long transitions period of the so called 
"preheating". The non-equilibrium character of 
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superheavy particle production implies strong de- 
pendence on the concrete physical processes that 
can take place at different periods of preheating 
stage. 

It was shown in |S] that the parametric reso- 
nance 6 in the end of inflation at t ~ 1/H en d, 
when preheating begins, can lead to intensive in- 
flaton field decay, in which superheavy particles 
with the mass m < 10H enc i can be produced. 
Here H en d ~ 10 13 GeV is the Hubble constant 
in the end of inflation. The calculations [5] of 
primordial concentration of such superheavy par- 
ticles exhibit strong dependence onm/H and cor- 
respond to a wide range of their modern densities 
up to fix ~ 0.3. 

Superheavy particles can be created in the end 
of preheating, when reheating takes place at t ~ 
1/H r (H r being the Hubble constant in the pe- 
riod of reheating), if the quanta of inflaton field 
contain these particles among the products of de- 
cay. The modern density of superheavy particles 
is then given by 



J-RD W0 
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where T r ~ (HrMp) 1 / 2 is the reheating temper- 
ature, Trd ~ lOeV is the temperature in the 
end of radiation dominance stage and in the be- 
ginning of the modern matter dominated stage, 
rrirf, > 2m is the mass of the inflaton field quan- 
tum and Br(X) is the branching ratio of super- 
heavy particles production in inflaton decay. The 
condition fix < 0.3 constrains the branching ra- 
tio as 



Br(X) < 0.1 



Trd m<t, 
T r m 



(3) 



If inflation ends by the first order phase tran- 
sition, bubble wall collisions in the course of true 
vacuum bubble nucleation can lead to formation 
of primordial black holes (PBH) with the mass 
M ~ Mpjj^ (see ^7] for review). Successive 
evaporation of such black holes at 
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is the source of superheavy particles, when the 
temperature of PBH evaporation, increasing with 



the loss of mass as Tpbh ~ Mp/M, reaches m. 
If ax is the fraction of PBH mass, evaporated 
in the form of considered superheavy particles, 
the relationship between the probability of PBH 
formation w and fix is given by 



fix = 



Ted 



a x w, 



(5) 



for dust-like (MD) expansion law at preheating 
stage and 



x 



(MpH en d) 
Trd 



1/2 



-axw, 



(0) 



for relativistic (RD) expansion at the stage of 
preheating. In the latter case corresponding to 
the Eq.@, the condition fix < 0.3 leads to 
w < 3 • 10~ 25 /ax- Creation of mini black holes 
with such a low probability does not imply first 
order phase transition after inflation, but it is pos- 
sible even from Gaussian "tails" (see ^2] for re- 
view) of nearly flat ultraviolet spectra, that are 
strongly disfavored but still not excluded within 
the uncertainty of the recent WMAP measure- 
ments of CMB anisotropy JH| • 

The presence of additional dynamically sub- 
dominant fields at the inflationary stage can 
strongly modify at the small scales the simple pic- 
ture of nearly flat power spectrum of density fluc- 
tuations. It also leads to the possibility of super- 
heavy particle production in the decay of quanta 
of such field, 4>, at the preheating stage. The re- 
lationship between fix and the relative contribu- 
tion r of the field, p^, into the total density pt D t, 
r = P4>/ptot in the period of decay, at r ~ l/Hd, is 
given by Eqs. Q-©, in which ax has the mean- 
ing of the branching ratio for superheavy particle 
production (multiplied by the factor ~ m/m^, 
m<f, is the mass of </>, in case of relativistic decay 
products) and H en d is substituted by Hd- If <f> de- 
cays due to gravitational interaction, Hd is equal 
to the probability of decay, T, given by Eq.QJ, 
H d = T ~ m^(m^/Afp) 4 . In general, for H d = T, 
the period of derelativization of the relativistic 
decay products with the energy e ~ 3> m 
corresponds to H ~ (m/m^) 2 ^. 
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3. Primordial pairing of superheavy parti- 
cles 

Note, first of all, following [5], that in quantum 
theory particle stability reflects the conservation 
law, which according to Noether's theorem is re- 
lated with the existence of a conserved charge, 
possessed by the considered particle. Charge con- 
servation implies that particle should be created 
together with its antiparticle. ft means that, be- 
ing stable, the considered superheavy particles 
should bear a conserved charge, and such charged 
particles should be created in pairs with their an- 
tiparticles at the stage of preheating. 

Being created in the local process the pair is 
localized within the cosmological horizon in the 
period of creation. If the momentum distribu- 
tion of created particles is peaked below p ~ mc, 
they don't spread beyond the proper region of 
their original localization, being in the period 
of creation I ~ c/H, where the Hubble con- 
stant H at the preheating stage is in the range 
H r < H < H en d- For relativistic pairs the re- 
gion of localization is determined by the size of 
cosmological horizon in the period of their derel- 
ativization. In the course of successive expansion 
the distance I between particles and antiparticles 
grows with the scale factor, so that after reheat- 
ing at the temperature T it is equal to 



1 H ' T 



(7) 



The averaged number density of superheavy 
particles n is constrained by the upper limit on 
their modern density. Say, if we take their max- 
imal possible contribution in the units of crit- 
ical density, fix, not to exceed 0.3, the mod- 
ern cosmological average number density should 



be 



10 



-20 W^GeV fix 



(being n = 4 • 



m 0.3 

1Q -22 W^GeV U2L T 3 j Q th j t h = C = k = 1 dt 
m 0.3 

the temperature T) . It correponds at the temper- 
ature T (at the redshift z) to the mean distance 
{l a ~ n" 1 / 3 ) equal to 



l s w 1.6 • 10 7 ( 



m U/3 f 0- 3 Nl/3 1 

10 14 GeV' [ n x ' T 



4 - 6 ' 106( il)^7 )1/3( S )1/3(1 + zrlcw - (8) 



One finds that superheavy nonrelativistic par- 
ticles, created just after the end of inflation, 
when H ~ H en d ~ 10 13 GeV, are separated 
from their antiparticles at distances more than 
4 orders of magnitude smaller, than the aver- 
age distance between these pairs. On the other 
hand, if the noncquilibrium processes of super- 
heavy particles creation (such as decay of infla- 
ton) take place in the end of preheating stage, 
and the reheating temperature is as low as it is 
constrained from the effects of gravitino decays 
on 6 Li abundance (T r < 4 • 10 6 GeV |17ll9p . the 
primordial separation of pairs, given by Eq(|7)>. 
can even exceed the value, given by Eq.JHJ). It 
means that the separation between particles and 
antiparticles can be determined in this case by 
their averaged density, if they were created at 
H < H s ~ 10~ 15 ■ M P (!2^Z) 2 / 3 (£f ) 2 / 3 ~ 

104(^)2/3(lol^el/)2/3 GeV /" 

If the considered charge is the source of a long 
range field, similar to the electromagnetic field, 
which can bind particle and antiparticle into the 
atom-like system, analogous to positronium, it 
may have important practical implications for 
UHECR problem. The annihilation timescale of 
such bound system can provide the rate of UHE 
particle sources, corresponding to UHECR data. 



4. Formation of primordial bound systems 
from primordial pairs 

The pair of particle and antiparticle with oppo- 
site gauge charges forms bound system, when in 
the course of expansion the absolute magnitude 
of potential energy of pair V = ^f- exceeds the ki- 

2 

netic energy of particle relative motion = ■ 
The mechanism is similar to the proposed in |2()| 
for binding of magnetic monopole-antimonopolc 
pairs. It is not a recombination one. The binding 
of two opposite charged particles is caused just 
by their Coulomb-like attraction, once it exceeds 
the kinetic energy of their relative motion. 

In case, plasma interactions do not heat super- 
heavy particles, created with relative momentum 
p < mc in the period, corresponding to Hubble 
constant H > H s , their initial separation, being 
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of the order of 

m = ( ^), ( 9) 

experiences only the effect of general expansion, 
proportional to the inverse first power of the scale 
factor, while the initial kinetic energy decreases as 
the square of the scale factor. Thus, the binding 
condition is fulfilled in the period, corresponding 
to the Hubble constant H c , determined by the 
equation 

(10) 



( H_ ]1/2 = 



P 



(—V - 
H r 2m 2 a„H' 



where H is the Hubble constant in the period 
of particle creation and a y is the "running con- 
stant" of the long range interaction, possessed by 
the superheavy particles. If the local process of 
pair creation does not involve nonzero orbital mo- 
mentum, due to the primordial pairing the bound 
system is formed in the state with zero orbital 
momentum. The size of bound system exhibits 
strong dependence on the initial momentum dis- 
tribution 
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P 



= 2- 



2a y m 3 H 2 to/3 2 



8-icr 6 (50a,)( lu ~ ev )(^—y cm , (11) 



rn 



P 



(12) 



where 

2a v mH 
(3 = — — 2 — . 

P 

In principle, it facilitates the possibility to fit 
UHECR data in the framework of hypothesis of 
bound system annihilation in the halo of our 
Galaxy. 

Indeed, the annihilation timescale of this bound 
system can be estimated from the annihilation 
rate, given by 

„2 



Qt 

w ann =| 1/(0) | 2 (av) ann ~ l- 3 -%C y , (13) 

to 



where the "Coulomb" factor C y arises similar to 
the case of a pair of electrically charged particle 
and antiparticle. For the relative velocity - < 1 
it is given by [21] 

C y = (14) 



Finally, taking v/c ~ /?, one obtains for the anni- 
hilation timescale 

1 , P NlO/ m N6 1 4 

~ (tf) — ~ — 



For H end > H > H s 
timescale equals 



(15) 



7TTO/3 5 

the annihilation 



io*(±°_^) ( i° )«., 

TO p 



(16) 



being for p ~ toc, a y = and to = 10 14 GeV in 
the range from l(T 26 s up to 10 19 (g|-) 10/3 s. The 
size of a bound system is given by 



^ = 8.10-(50 % )( 10l4GeF - 10 " 2 ^ 



rn 



-)(- 



fcm, (17) 



ranging for 2 • 10 _ 10 < §f < 1 from 7 • 10~ 7 cm 
to 6 • 10 _3 cm. One can obtain from Eqs. (|16fl - 
l|17|l the approximate relationship between T ann 
and l c , given by 



(- 



lc. 



5/2 



10 1( V v 10- 7 cto/ ' ^ 

Provided that the primordial abundance of su- 
perheavy particles, created on preheating stage 
corresponds to the appropriate modern density 
fix < 0.3, and the annihilation timescale exceeds 
the age of the Universe tjj = 4 • 10 17 s, owing to 
strong dependence on the parameter (3, the mag- 
nitude 

rx = 03 ~ x (19) 

can easily take the value rx = 2- 10~ 10 , which was 
found in to fit the UHECR data by superheavy 
particle decays in the halo of our Galaxy. It takes 
place, provided that 

If the effective production of superheavy parti- 
cles takes place in the end of preheating stage at 
H < H s , their initial separation is determined by 
the min{l(H), l s }, where 1(H) is given by Eq. 
and l s is determined by their mean number den- 
sity (compare with Eq.©) 



Z s « 3 • 10 7 (- 



10 14 GeV 



\l/3/_°_£\l/3/ 1 



M P H 



1/2 



6 



10" 4 ( 



10 14 GeV 



^ ,/3 (g) 1/3 (^^) 1/2 -.(21) 



In the case of late particle production (i.e. at 
H < H s ) the binding condition can retain the 
form 1|1U[I. if 1(H) < l s . Then the previous esti- 
mations <|11|) - H15fl are va lid- 

5. Formation of primordial bound systems 
without primordial pairing 

In the opposite case of late particle production, 
when 1(H) > l s , the primordial pairing is lost and 
even being produced with zero orbital momen- 
tum particles and antiparticles, originated from 
different pairs, in general, form bound systems 
with nonzero orbital momentum. The size of the 
bound system is in this case obtained from the 
binding condition for the initial separation, de- 
termined by Eq.J2U, and it is equal to 



10 



15 



2a v M P y 10 14 GeV 



) 2/3 (^) 2/:i ^^77' 



JL 

mc 



H' 



2-icr 6 ( 



_j|_ )2/ . ( o£ )2/ . ( ^! )cm . (22: 



The orbital momentum of this bound system 
can be estimated as M ~ mvl c and the life- 
time of such bound system is determined by the 
timescale of the loss of this orbital momentum. 
This timescale can be reasonably estimated with 
the use of the well known results of classical prob- 
lem of the falling down the center due to radiation 
in the bound system of opposite electric charges 
ei and e 2 with masses m\ and 1112, initial orbital 
momentum M and absolute value of the initial 
binding energy E (see e.g. [22] 1 

„3 ]\/r5 



c 6 M b ei e 2 a 
f a y (2Efi 3 y/ 2[ mi m 2 ' ' 



(23) 



((ml) 1/2 + (2M 2 Ef' 2 ) 

Here a = mi | T " 2 is the reduced mass. Putting 
into Eq.|J2JJ| M = fj,vl c , E = \iv 2 j2, and with the 
account for 2M 2 E ~ fia 2 one obtains the lifetime 
of the bound system as 



r = 



64-7T a 2 



4-10 2U ( 



L 



-) 3 ( 



10" 6 cm' y 10 14 GeV 



) 2 (50a y )- 2 yr.(24) 



Using the Ea. H13|l and the condition 1(H) > l s , 
one obtains for this case the following restriction 



rx 



fix tu 
"03^ 



<3.10-»&«(^K)».(25) 

0.3 771 



Note that the condition (|25|l admits fix — 0.3, 
when superheavy particles dominate in the dark 
matter of the modern Universe. 

The gauge U(l) nature of the charge, possessed 
by superheavy particles, assumes the existence of 
massless U(l) gauge bosons (y-photons) mediat- 
ing this interaction. Since the considered super- 
heavy particles are the lightest particles bearing 
this charge, and they are not in thermodynami- 
cal equilibrium, one can expect that there should 
be no thermal background of y-photons and that 
their non equilibrium fluxes can not heat signifi- 
cantly the superheavy particles. 

6. Primordial bound systems in the case of 
plasma heating 

The situation changes drastically, if the super- 
heavy particles possess not only new U(l) charge 
but also some ordinary (weak, strong or elec- 
tric) charge [SJ- Due to this charge superheavy 
particles interact with the equilibrium relativistic 
plasma (with the number density 77 ~ T 3 ) and 
for the mass of particles 777 < a 2 Mp the rate of 
heating 



navAE ~ a 2 — 
777 



(26) 



is sufficiently high to bring the particles into ther- 
mal equilibrium with this plasma. Here a is the 
running constant of the considered (weak, strong 
or electromagnetic) interaction. 

Plasma heating causes the thermal motion of 
superheavy particles. At T < rn( a2 ™ Ip ) 2 their 
mean free path relative to scattering with plasma 
exceeds the free thermal motion path, so it is 
not diffusion, but free motion with thermal ve- 
locity vt that leads to complete loss of initial 
pairing, since vxt formally exceeds l s at T < 

10-10M P (|f) 2 / 3 (±2±l^)5/3_ 
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In the case, the interaction with plasma keeps 
superheavy particles in thermal equilibrium, po- 
tential energy of charge interaction V — 

is less, than thermal energy T for any a y < 
3 • 107(0^)1/3( ^^ )1/3^ So binding condition 

V > Tki n can not take place, when plasma heat- 
ing of superheavy particles is effective. 

For electrically charged particles it is the case 
until electron positron pairs annihilate at T e ~ 
lOOkeV (see [2D] ) an d for colored particles until 
QCD phase transition at Tqcd ~ 300MeV. In the 
latter case colored superheavy particles form su- 
perheavy stable hadrons, possessing U(l) charge. 
For weakly interacting particles after electroweak 
phase transition, when Eq. (|20[) is not valid, neu- 
trino heating, given by navAE ~ G 2 F — , is suffi- 
ciently effective until T w w 20GeV. AkT < T N , 
where N = e,QCD,w respectively, the plasma 
heating is suppressed and superheavy particles go 
out of thermal equilibrium. 

In the course of successive expansion kinetic 
energy of superheavy particles falls down with the 
scale factor a as a a 2 , and the binding condition 
is reached at T c , given by 

^ = ^3.10^(^3(1^^)1/3. ( 27) 

0.3 m 

For electrically charged particles, forming after 
recombination atom-like states with protons and 
electrons, but still experiencing the Coulomb- like 
attraction due to non-compensated U(l) charges, 
the binding in fact does not take place to the 
present time, since one gets from Eq.l|27|l T c < 
IK. Bound systems of hadronic and weakly in- 
teracting superheavy particles can form, respec- 
tively, at T c ~ 0.3eV and T c « 20eV. 

The size of the bound system is then given by 

lc = 10 15 (^) 2 / 3 (— ^— ) 2 / 3 ^- » (28) 
\lx 10 14 GeV / a.y±N 

what even for weakly interacting particles ap- 
proaches a half of meter (30 m for hadronic 
particles!). It leads to extremely long annihi- 
lation timescale of these bound systems, that 
can not fit UHECR data. Moreover, being 



extremely weakly bound, they should be dis- 
rupted almost completely, colliding in Galaxy. 
So, for bound systems of weakly interacting su- 
perheavy particles, navtjj ~ 10 14 where n = 
3 ■ 10~ 15 cm~ 3 wl4 £ eV §f is the number density of 
bound systems, a ~ irl 2 and their relative velocity 
v ~ 3T0 7 cm/s. It makes impossible to realise the 
considered mechanism of UHECR origin, if the 
superheavy U(l) charged particles share ordinary 
weak, strong or electromagnetic interactions. 

7. Evolution of bound systems in the 
Galaxy 

Superheavy particles, as any other form of 
CDM should participate gravitational clustering 
and concentrate to the center in the course of 
Galaxy formation. There are several factors in- 
fluencing the evolution of bound systems in the 
Galaxy. 

If the size of primordial bound systems l c > 
3 • 10~ 6 cm(^|- w iF GeV ) 1/2 their collision rate in 
the vicinity of Solar system navtjj > 1. It can 
take place, provided that 

(50a,) 4 (H^)n/ 5 (^| )9 / 5 > 6 . i -« (29) 

Since the binding energy of bound systems 
(E b < 5 • 10- 27 mc 2 (^- 10l4 w f y ) 2 / 5 ) is much less 
than the kinetic energy of their motion in the 
Galaxy (T = mv 2 /2), the bound systems should 
be disrupted in such collisions. On the other 
hand, large collision cross section a ~ irl 2 corre- 
sponds to the momentum transfer of the order of 
Ap ~ m[3c. Such small momentum transfer leads 
both to disruption of the bound system on free 
superheavy particles and with the same order of 
probability to the reduction of their size down to 
I ~ Z c /4. If the bound systems are within gravita- 
tionally bound cluster (Galaxy, globular cluster, 
CDM small scale cluster), both the free particles 
and contracted bound systems with the size I re- 
main in it and the collisions both between bound 
systems and between free particles and bound 
systems continue with the cross section a ~ irl 2 . 

In vicinity of massive objects with the mass M 
tidal effects lead to disruption of bound systems 
with the size I at the distance r, corresponding 
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to tidal force energy 



GMm I 



-, exceeding the 



Galaxy, is given by 



binding energy ~ a y /l. So, in the vicinity of 
a star with the mass M — M@ the bound sys- 
tems are disrupted at the distances, smaller than 
r <i = \ G t Im l- Tidal effects disrupt bound sys- 

V v 

terns in the regions of enhanced stellar density, 
um, if TiMKr^vtu > 1. Taking for the center of 
Galaxy um ~ 10 6 Mq/pc 3 , one finds that bound 
systems with the size I > 5 • 10~ 6 cm should be 
disrupted there due to tidal effects. 

As it was recently shown in 123], tidal effects 
strongly influence the formation and mass distri- 
bution of small scale CDM clusters, what should 
also take place for small scale clusters of pri- 
mordial bound systems in the Galaxy. On the 
other hand, clustering of primordial bound sys- 
tems may play important role in the explanation 
of observed clustering of UHECR events in the 
framework of the proposed mechanism. It may 
be easily estimated that if bound systems are 
clustered around globular cluster, their disrup- 
tion due to stellar tidal effects is negligible. 



8. Space distribution and clustering of 
UHECR events 

For the initial number density and size of bound 
systems, corresponding to navtu > 1, most of 
bound systems disrupt on the free particles, but 
sufficiently large fraction of them ~ (3 C /Pu ~ 
(lu/lc) 1 ^ 2 acquires the size lu, at which navtu ~ 
1. The relative amount of bound systems with 
smaller size I < Ijj is of the order of ~ (l/ljj) 2 , 
if their annihilation timescale t > tu and of the 
order of ~ (l/lu) 2 jr; for r < tjj . Annihilation 
of superheavy particles in bound systems with 
the smaller size is more rapid, what increases 
the production rate of such UHECR source as 
compared with the case of superheavy decaying 
particles with the fixed lifetime. This effect of 
the self- adjustment of bound systems annihila- 
tion leads to the peculiar space distribution of 
UHECR sources, corresponding to this mecha- 
nism. 

The decay rate density, q, of metastable parti- 
cles with lifetime r for the number density n(R), 
depending on the distance R from the center of 



n(R), 



(30) 



Owing to self adjustment of bound systems anni- 
hilation the density of UHECR production rate 
for the initial number density n(R) of primordial 
bound systems with initial annihilation timescale 
r has the order of the magnitude 



q 



, n(R) t 3/5 r . 1/s ^ 

T Tu tu 



(31) 



Since tu oc l^/ 2 oc n(i?) -5 / 4 the self-adjustment 
of bound systems leads to stronger radial depen- 
dence q oc n(R) 7 / 4 , thus sharpening the con- 
centration of UHECR sources to the center of 
Galaxy. 

In the case of late particle production with 
1(H) > l s , considered in the Section 5, the de- 
pendence I oc /3 _1 , given by Ea. (|22|l . and t oc 
I 3 oc f3~ 3 , given by Ea. l|24Jl . results in the change 
of Eq.JHU) by 

?=(^)(-) 2/3 (f) 1/3 . 

T TU tu 

Since tu oc l v oc n(R)~ 3 / 2 in this case, the radial 
dependence q oc n(R) 3 / 2 provides the principal 
possibility to distinguish this case from the case 
1(H) < l s . 

It was noticed in [23 that clustering of UHECR 
events, observed in AGASA experiment |25|. 
can be explained in the model of superheavy 
metastable particles, if such particles with the 
mass m ~ 10 GeV form clusters in the Galaxy 
with the mass M ~ 5 • Mm tfStt — ■ 

For the cluster of N metastable particles with 
lifetime r the decay rate, P, is given by 



p = (-). 

T 



(32) 



Owing to self adjustment of bound systems anni- 
hilation the UHECR production rate for the clus- 
ter of N primordial bound systems with initial 
annihilation timescale r reaches the order of the 
magnitude 



^=(-)(-) 3/5 (f) 1/5 - 

T TU tu 



(33) 
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The enhancement of UHECR production rate 
due to self-adjustment of bound systems facil- 
itates the possibility to explain clustering of 
UHECR events in the proposed mechanism. Say, 
instead of cluster with the mass of ~ 5 • 10 6 M Q 
of particles with the mass m ~ 10 14 GeV and the 
lifetime r ~ 10 16 yr, corresponding to the consid- 
ered in possible explanation for the cluster- 
ing of UHECR events, observed by AGASA, it 
is sufficient to have the mass of ~ 3 • 10 3 M Q in 
the cluster of primordial bounds systems with the 
same initial annihilation timescale r « 10 16 yr. 

In the case of late particle production with 
1(H) > l S) for p ~ me, the condition H < H s 
constrains the annihilation timescale 



m 



10 14 GeU 



) 9 (50a y )" 



= 3 - 1 ° 1 V( I ^) 9 (50%)- 3 , 

the size of the bound system I being related with 
r by Ea. H24|) . In this case, owing to the effect of 
self- adjustment, clustering of UHECR events can 
be also reproduced, e.g. in cluster with the mass 
M ~ 10 5 M Q and number density n ~ I0 _10 cm~ 3 
of superheavy particles with mass m ~ 10 13 GeU 
and initial annihilation timescale of their bound 
systems r = 4 • 10 ls yr (for metastable particles 
with the same mass and lifetime the cluster with 
the mass 4 ■ 10 8 M Q is needed). Taking into ac- 
count the possibility of dominance of superheavy 
particles in the modern CDM, mentioned above 
for the considered case, the estimated parameters 
of such cluster seem to be reasonable. 

Owing to the effect of self-adjustment the in- 
tial annihilation time decreases in the dense re- 
gions. On the other hand, disruption of bound 
systems in collisions leads to the decrease of their 
actual amount in the modern Universe as com- 
pared with their primordial abundance. It leads 
to the corresponding corrections in the conditions 
(|2Uj) and J25J|. Provided that the inequality 
is valid and the collisions of bound systems are 
significant, the annihilation timescale r should be 
corrected by the effect of self-adjustment and one 
should substitute fix by flt s < flj;, where Qb s 
is the averaged concentration of bound systems, 
surviving after collisions. 



9. Annihilation into ordinary particles 

To be the source of UHECR the products of 
superheavy particles annihilation should contain 
significant amount of ordinary particles. On the 
other hand, it was shown in the present paper 
that to be the viable source of UHECR the 
considered particles should not possess ordinary 
strong, weak and electromagnetic interactions. 
Their interaction with ordinary particles, giving 
rise to UHECR production in their annihilation 
should be related to the superhighenergy sec- 
tor of particle theory and/or physics of inflation 
and preheating. The selfconsistent treatment of 
this problem should involve the realistic particle 
physics model, reproducing the desired features 
of inflationary scenario and giving detailed pre- 
dictions for physical properties of U(l) charged 
superheavy particles. It may be expected that 
in such models, there can exist superheavy bo- 
son (Y), interacting both with superheavy and 
ordinary particles. If it's mass is of the order of 
my > to the annihilation channel into ordinary 

2 

particles with the cross section a ~ -^-j-rn 2 will 

m w 

be of the order of the cross section of the two 
y-photon annihilation channel. 

With the account for the invisible yy mode of 
annihilation, as well as taking into account effects 
of self-adjustment and destructions of bound sys- 
tems in collisions, one has to re-define the mag- 
nitude rx, given above by Eq. I|19|) . as 



(34) 



0.3 r eff 

where B a is the branching ratio of annihilation 
channels to ordinary particles and Qbs < &x is 
the modern cosmological density of bound sys- 
tems. In the Eq. $Q r e ff = (t x T^iy) 1 / 5 and 
T eff — (TxTutu) 1 ^ 3 for the cases of early and late 
particle production, respectively. 

Note that at B ~ I CP 5 , when annihilation 
to ordinary particles is strongly suppresssed, the 
case fix ~ 0.3 is possible also for early particle 
production. 

Neutrino channel may strongly dominate in 
the annihilation to ordinary particles, so that 
B Q rs B u and the channels to other ordinary 
particles are strongly suppressed (B op -C I0~ 5 , 
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where op = q,l,j,g,W, Z,h). Then annihilation 
of bound systems can not be direct local (galactic) 
source of UHECR, but it can provide the source 
of UHE neutrinos for the Z-Shower mechanism of 
UHECR origin. 

10. CONCLUSIONS 

The combination of the constrains on the con- 
ditions of particle creation in the early Universe 
and on the effective production of UHECR puts 
additional constrains on the parameters of the 
proposed mechanism, which we plan to consider 
in the framework of specific models of particle 
theory, underlying the scenarios of very early Uni- 
verse. The evolution of primordial bound systems 
should be also analyzed on the base of such mod- 
els. One, however, can make the general conclu- 
sion that the two principal types of bound sys- 
tems are possible, originated from (i) "Early par- 
ticle production" , when primordial pairing essen- 
tially determines the formation of bound systems 
and from (ii) " Late particle production" , when 
the primordial pairing is not essential for bound 
system formation. 

In the both cases bound systems can dominate 
in the modern CDM, but the conditions for such 
dominance are different. In the case (i) anni- 
hilation of bound systems with f2{, s ~ 0.3 can 
reproduce UHECR events, if the branching ra- 
tio for annihilation to ordinary particles is small 
(B a ~ 10 -5 ), whereas in case (ii) this branching 
ratio should maximally approach to 1. 

The possibility of bound system disruption in 
their collisions in galaxies is specific for the con- 
sidered mechanism, making it different from the 
models of decaying and annihilating superheavy 
particles. If effective, such disruption results in 
a nontrivial situation, when superheavy particles 
can dominate in the modern CDM, while the 
UHECR sources represent a sparse subdominant 
component of bound systems, surviving after dis- 
ruption. 

Self-adjustment of bound systems annihilation 
in the Galaxy sharpens their concentration to the 
center of Galaxy and increases their UHECR pro- 
duction rate in clusters. It provides their differ- 
ence from the case of metastable particles. This 



property, however, crucially depends on the prob- 
ability of contraction of bound systems in the 
course of collisions. More detailed analysis of such 
collisions is needed to prove this result. If proven 
to be really specific for the considered type of 
UHECR sources, it will be principally possible to 
distinguish them from other possible mechanisms 
|lf)lftlS6| in the future AUGER and EUSO exper- 
iments. 

Pair correlation, considered in the present note, 
takes place, if the local process of superheavy par- 
ticle creation preserves charge conservation. This 
condition has serious grounds in the case of a local 
U(l) gauge charge, similar to electric charge, but 
it may not be the case for global charge, say, for 
mechanisms of R-parity nonconservation due to 
quantum gravity wormhole effects ■ The crucial 
physical condition for the formation of primordial 
bound systems of superheavy particles is the ex- 
istence of new strictly conserved local U(l) gauge 
symmetry, ascribed to the hidden sector of par- 
ticle theory. Such symmetry can arise in the ex- 
tended variants of GUT models (see e.g. ^Jj for 
review), in heterotic string phenomenology (see 
|21j and references wherein) and in D-brane phe- 
nomenology ■ Note, that in such models the 
strictly conserved symmetry of hidden sector can 
be also SU(2), what leads to a nontrivial mech- 
anism of primordial binding of superheavy parti- 
cles due to macroscopic size SU(2) confinement, 
as it was the case for "thetons" |28j . 

The proposed mechanism is deeply involved 
into the details of the hidden sector of particle 
theory, what may seem rather artificial and fine- 
tuned. However, the necessary combination of 
conditions (superheavy stable particles, possess- 
ing new strictly conserved U(l) charge, existence 
of their superheavy Y-boson interaction with or- 
dinary particles, nontrivial physics of inflation 
and preheating) can be rather naturally realised 
in the hidden sector of particle theory. In this 
aspect, the proposed mechanism offers the link 
between the observed UHECRs and the predic- 
tions of particle theory, which can not be tested 
by any other means and on which the analysis 
of primordial pairing and binding can put severe 
constrains. 

Even so, while we may agree that the number 
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and sequence of the assumption of present sce- 
nario may sound artificial and ad hoc, (maybe at 
the same level of topological decay lifetime) we 
have taken into account a large number of astro- 
physical and cosmological bounds narrowing the 
parameter window into a very severe and frag- 
ile regime which may soon survive (or not) fu- 
ture theoretical self-consistence and experimen- 
tal test. Indeed, if HIRES and AGASA data will 
converge to a GZK cut off with no spectra exten- 
sions to Grand Unified energies or, in a different 
scenario, in case of more evidence for UHECR 
clustering to BL Lac sources compatible only to 
Z-Showcr model |l()lllj . the model may be con- 
sidered as an untenable direct local solution of 
UHECR puzzle. However, for suppressed chan- 
nels of annihilation to quarks, charged leptons, 
gauge and Higgs bosons, annihilation of super- 
heavy particles in bound systems can provide the 
effective sourse of UHE neutrinos, thus playing 
important role in the UHECR production by Z- 
Shower mechanism. It should be remind that 
UHE neutrinos of all flavours will be produces 
by such heavy particle annihilation leading to im- 
portant signals in new generation UHE neutrino 
telescopes based on Horizontal Tau Showering (or 
Earth Skimming Neutrinos) (2H| (see review and 
some answers to the recent arguments (SO) against 
Z-Shower mechanism in |31|L 

If viable, the considered mechanism makes 
UHECR the unique source of detailed informa- 
tion on the possible properties of the hidden sec- 
tor of particle theory and on the physics of very 
early Universe 
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